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We develop the "triangulated" version of loop quantum cosmology, recently introduced in the liter- 
ature. We focus on the "dipole" cosmology, where space is a three-sphere and the triangulation is 
formed by two tetrahedra. We show that the discrete fiducial connection has a simple and appealing 
geometrical interpretation and we correct the ansatz on the relation between the model variables 
and the Friedmann-Robertson- Walker scale factor. The modified ansatz leads to the convergence 
of the Hamiltonian constraint to the continuum one. We then ask which degrees of freedom are 
captured by this model. We show that the model is rich enough to describe the (anisotropic) Bianchi 
IX Universe, and give the explicit relation between the Bianchi IX variables and the variables of 
the model. We discuss the possibility of using this path in order to define the quantization of the 
Bianchi IX Universe. The model contains more degrees of freedom than Bianchi IX, and therefore 
captures some inhomogeneous degrees of freedom as well. Inhomogeneous degrees of freedom can 
be expanded in representations of the SU(2) Bianchi IX isometry group, and the dipole model cap- 
tures the lowest integer representation of these, connected to hyper-spherical harmonic of angular 
D ■ momentum j = 1. 
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D" 1 I. INTRODUCTION 

Loop quantum cosmology (LQC) [l[ is the most remarkable application of the loop approach to quantum gravity. Its 
main result is a robust indication that the cosmological singularity that appears in classical general relativity is removed 
£Nj . by quantum effects [2J . While the physical basis of this results is explicitly grounded on the physical discreteness of 
quantum geometry which is predicted Q by full loop quantum gravity (LQG) Q, the precise relation between LQC 
and LQG, on the other hand, has not been fully clarified yet To shed light on this issue, a finite dimensional 
truncation of LQG, which leads naturally to LQC in the Born-Oppenheimer approximation, was introduced in 
The idea is to fix a coarse triangulation of physical space and consider a discretization and a quantization of general 
relativity on this triangulation. This procedure leads to a truncated version of LQG, interpreted as a description of 
a finite number of large-scale (namely cosmological) degrees of freedom. In the simplest version of the theory, called 
the "dipole" cosmology, compact physical space is triangulated with two tetrahedra. It was shown in [6| that this 
model leads to a LQC-like dynamics in the Born-Oppenheimer approximation, where the scale factor plays the role 
of the "heavy" degree of freedom. The discrete dynamics of LQC is recovered in this manner without recurring to 
the "area gap" argument. 

The main purpose of this work is to analyze this dipole cosmology more in detail, and in particular, interpret the 
geometrical meaning of its other ("light") degrees of freedom. It has been suggested that these could just capture 
some space anisotropies. We show that the model captures indeed all the anisotropic degrees of freedom, that is, the 
degrees of freedom of Bianchi IX. This may be of particular interest since Bianchi IX describes a generic space-time 
near the cosmological singularity 7] . We write an explicit relation between the Bianchi IX degrees of freedom and 
those of the model. We discuss the possibility of using this relation in order to construct a loop quantization of the 
Bianchi IX cosmological model. 

However, the model captures a richer dynamics than Bianchi IX: it includes inhomogeneous degrees of freedom 
as well, and therefore truly represents a "first step towards inhomogeneity" in LQC. These degrees of freedom can 
be identified as the lowest order term of a tensor harmonic expansion, in terms of Wigner functions, of the spatial 
geometry (see Ref. [8]). The physical relevance of this expansion is of interest not only in cosmology, but also in 
relation with the study of the LQG n-points correlation functions, which is based on a similar approximation Hoj. 
In the course of our analysis, we obtain also two other results. First, we show that the discrete version of the 
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fiducial connection used to define the theory, which was defined in a rather ad-hoc manner in [6J, has in fact a neat 
and appealing geometrical interpretation on the triangulated space. Second, following we correct the naive Bohr- 
Oppenheimer ansatz used in Q to relate the model variables to the scale factor. We show that using a different ansatz 
the Friedman equation can be obtained explicitly in the appropriate limit, from the discrete Hamiltonian constraint. 
Therefore the model agrees with the standard formulation of LQC both at the classical and quantum levels. The 
triangulated Hamiltonian constraints can be interpreted as a different version of the LQC constraint. 

We recall that although most of the LQC literature focuses on the isotropic sector [l^, EH, homogeneous [l4T - fl(| 
and inhomogeneous [171 . Il8j models have been studied as well. 

The paper is organized as follows. In Section QI] we review the geometry of the Bianchi IX Universe. In Section ITTT1 
we review the dipole cosmology model. In Section HVl we relate the the dipole cosmology with Bianchi IX. In Section 
IVl we study the classical limit of the model and show that it reproduces the classical theory at small curvature. In 
Section [VTl we study the way anisotropies are described in the model. Finally, inhomogeneous degrees of freedom are 
introduced in Section Ivm 

We set 8ttG/3 = 1 = c, where G is the Newton gravitational constant and c the speed of light. We use small 
Latin letters a, b, c, ... as spatial indices on the three-dimensional surfaces, and capital Latin letters /, J, K... to label 
internal indices of su(2)-algebra-elements, which are raised and lowered with the identity matrix. Sum over internal or 
external indices is meant whenever those are repeated, even if they are not paired, and do not appear inside brackets. 
When indices are within brackets, summation is not be intended. That is ajb 1 c 1 = aib 1 c 1 , but a/nb^'c^ = d 1 . 
We consider Euclidean gravity and we set for simplicity the Barbero-Immirzi parameter as ft = 1. 



II. BIANCHI IX IN A NUTSHELL 



In this Section we recall the basic properties of the Bianchi IX cosmological model, which plays a central role in 
the following. For a detailed discussion about the topology of this model we refer to the Appendix [Al 

The Bianchi IX model is the most general homogeneous model for a spatially compact Universe. It is also the basis 
of the classical description of a generic inhomogeneous space-time near the singularity, via the Belinski-Khalatnikov- 
Lifshitz (BKL) scenario 0]: as the cosmological singularity is approached spatial points (causal horizons) decouple 
dynamically and each of them evolves essentially independently as a Bianchi IX model [?J • Thus classical physics is 
well described in terms of this model near the cosmological singularity. 

Consider a four-dimensional metric space-time M. and let <f> : M. — > K. £3 E t be a diffeomorphism, where S t are 
the Cauchy surfaces foliating space-time. This is defined to be spatially homogeneous if for any i £ I and any two 
points p, 3 6 S( there exists an isometry of the space-time metric which takes p into q. A Bianchi IX cosmology is a 
homogeneous cosmology where E t has the topology of the three sphere S* 3 . 

In order to write homogeneous fields on S , it is convenient to have a homogeneous reference triad field. This can 
be constructed by exploiting the fact that S 3 can be identified with the group manifold of SU{2). This manifold 
carries the natural Cartan (flat) connection 

uj = g~ l dg = lj 1 ' tj = uj^Tidx a , (1) 

where g £ SU(2) and x a are three arbitrary coordinates on SU(2). Here 77 = o//(2z) is a basis in the su(2) Lie 
algebra and o~i are the Pauli matrices. This connection is left-invariant and satisfies the Maurer-Cartan structure 
equation 

dui 1 — j K ui' J A oj k = , (2) 

where tijK is the completely antisymmetric tensor. Denote ej = ejd a the corresponding dual vector field, with values 
in su(2), such that ejw a 7 = 8j . The Lie brackets of these vector fields read 

[ej, ej] = -e K jje K . (3) 

By identifying su(2) with an internal space, we can take ej and lu 1 as the definition of a triad and co-triad in space, 
which are invariant under the (left) action of SU(2) on itself. These fields can be used as "fiducial" co-frames and 
frames, as they naturally carry information about the homogeneity symmetry group. 
An explicit parametrisation for the co-frames can be given in terms of Euler angles as 

lui = cos ip d6 + sin ip sin 6 d<f) , (4) 
0J2 = sin ip dO — cos ip sin 9 dtp , 
W3 = dip + cos 9 d<f> . 
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Here 9 and <\> are in the range 9 £ [0,7r), <\> £ [0,27r), while for the Euler angle ip one must set ijj £ [0,4tt) in order 
to achieve a simply connected covering space with the topology of S 3 . This parametrization allows us to make easily 
contact with the definition of Cartan connections in terms of elements g € SU(2). These latter can in fact written as 

.9 = ~J ) , with 

a = eW+W cos(6»/2), = e -*W-*> sin(0/2) , (5) 

from which the above expression for the co- frames follows by means of |T]). One can also easily find the expression 
for the left invariant frame e" in this chart, as well as the right invariant vector fields (the Killing vector fields) which 
Lie drag the above frame and co-frame. 

Any homogeneous Ricmannian metric q a \, on the Cauchy surfaces can be then written in terms of its triadic 
projection qjj on internal space as 

q a b(x,t) = qu{t)u) I a {x)u)({x), (6) 

where qu is a 3 x 3 matrix, constant in space. In the isotropic case qu is a multiple of the identity, and it is 
proportional to the Killing-Cartan metric on SU(2). The space-time line-element can then be written in the form 
(taking the lapse function N = 1) 

ds 2 = -dt 2 + q u uj 1 ® u J . (7) 
The vacuum Einstein equations allow us to write qj j (t) in diagonal form (l9l I20I ] 

q ab = aj u&l , (8) 

where aj — ai(t) are the three (independent) scale factors which describe the anisotropy of the space slices £4. Thus 
homogeneity reduces the phase-space of general relativity to six dimensions. The inverse metric reads q ab = aj 2 ef ej 
and the spatial volume V is given in terms of the scale factors by V = 27r 2 ai 02 03. The closed Friedmann- Robertson- 
Walker (FRW) cosmological model is recovered as the particular case a% = 0,2 = 03. 

The Ashtekar-Barbero variables for the Bianchi IX cosmological model are hence written in terms of the Maurer- 
Cartan connection ([2]) and the dual triad, and are parametrized by three time- functions c 1 (t) and their duals p T {t) as 
(see for instance (ill. Hq) 

Ai = c (J)W W, Ef=pW U q x) , (9) 

in which enters the determinant of the co-triad u> = det(w^), namely the square root of the determinant of the metric 
of the three sphere. The connections cj and momenta p 1 parametrize the six-dimensional phase space, the symplectic 
two- form being = dci A dp 1 . The momenta are related to the metric variables by the relations 

Pi = |a 2 a 3 |sgn(ai), p 2 = \aia 3 \sgn(a 2 ), P3 = |aia 2 |sgn(a 3 ) , (10) 

and the connections are given in terms of triadic projection of Christoffel symbol components Ti and extrinsic curvature 
components Ki — — a//2 as c/ = Tj — K%. The Ti are here given in terms of the scale factors a/ as 

1 faj_ ok _ aj \ = l(JL,PL_ P J P K \ , m 
1 2\a K aj aja K ) 2 U J P A ' (p 1 ) 2 J ' 1 ' 

In the closed FRW model [l3[ (ai = a 2 = 03) the triadic projection of the Christoffel symbols becomes the constant 
r = 1/2 and one recovers the isotropic connection c — (d + l)/2 as well as the momentum \p\ = a 2 . Thus the 
cosmological singularity in Bianchi IX appears whenever aj = for some I. 



III. A MODEL FOR MERGING LQC IN LQG 



We briefly recall the construction of the model @ . 
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A. Classical theory 

Fix an oriented triangulation A n of the topological three-sphere, formed by n tetrahedra t glued by their triangles. 
Label the triangles with an index / ("/" for face) that runs from 1 to 2n (the number of faces is twice the number 
of tetrahedra). A group element Uf € SU(2) and a su(2) algebra element Ef are associated to each oriented triangle 
/. Given the face /~ x obtained inverting the orientation of /, we take the convection that its associated group and 
algebra elements read 

U f -i=Uj l , E f - 1 =-Uy 1 E f U f . (12) 

We use the notation Ef — Ejri and take Uf and Ef as the phase space variables of a dynamical system. The phase 
space of this model is that of a canonical lattice SU(2) Yang-Mills theory, i.e. the fundamental Poisson brackets are 
given by 

{Uf,U f ,}=0, {E I f,U r } = 6ff, T J U f , {E I f,E J f,}=-5 ff , e IJK Ef. (13) 

In other words, the phase space is the cotangent bundle of SU(2) 2n with its natural symplectic structure. The 
dynamics of the system is defined by two sets of constraints. The Gaufi (gauge) constraint (three constraints per 
tetrahedron) 

G t = (14) 
fet 

where the sum is over the four faces of the tetrahedron, and the Hamiltonian constraint 

C t ee Vf l J2 Tr\Uf f ,Ef,E f ] « . (15) 
ff'et 

In (fl"5|) the sum is over the couples of distinct faces at each tetrahedron, Uff = U fUJ, 1 and V 2 = Ti[EfEfEf»} because 
of (fl4)) . Here Vt can be interpreted as (proportional to) the volume of the tetrahedron t — thus we call V = J2t 
the total volume of space. The Hamiltonian constraint gives the good classical limit once the gauge constraint has 
been taken into account. In fact, at the first order the holonomy is U ~ exp J A ~ 1 — \"f\ 2 F + 0(|7| 4 A 2 ), where F 
is the curvature of the connection A and 7 denotes a loop. By means of the holonomy expansion, equation (|15p can 
be formally recast as follows 

V t C t = J2 ^[E r Ef] - | 7 | 2 Y, ^l F ff E f' E f\ ~ °- ( 16 ) 
ff'et ffet 

The former term in relation (jTHJ) vanishes because of the Gaufi constraint, i.e. ^2fj, £t Ti[EfEf] = 
^ r [(S/et E f)Cl2f'et Eft)] ~ 0, and thus the second term undergoes the expected continuum limit. We stress that 
this happens not only for small values of the length of the loop |7|, but it does for large values of | 7 | too, provided 
that \"/\ 2 F is small. 

This model can be regarded as a lattice approximation of the geometrodynamics of a closed Universe. To see 
this, consider real Ashtekar connection A ! a and electric field Ef, with their standard Poisson algebra (see for instance 
[H), on a three-dimensional surface S with the S 3 topology. Let A„ be a triangulation of S and A* a dual of the 
triangulation. Then, Uf is the parallel transport of the Ashtekar-Barbero connection along the link e/ of A* dual 
to the triangle /. In the same way, Ef is the flux $/ of the conjugate electric field across the triangle / (parallel 
transported to the center of the tetrahedron). The Poisson brackets of Uf and Ef are thus exactly the ones in (|13p . 
For instance, the third equality in (fT3")) is motivated by the parallel transport of to the center of each tetrahedron. 
Of course these variables transform via the gauge constraint (TT4")) . 

The constraint (1151) corresponds to the non-graph-changing version of the Hamiltonian constraint. Its quantization 
is an operator acting on an underlying spin-network state representing the graph dual to the triangulation taken in 
consideration on the spatial manifold. Whenever the triangulation is sufficiently fine, constraint (1151) represents an 
approximation the to Euclidean part of the Hamiltonian constraint Ti[F a i,E a E b ) / Vdct E ps 0. 

B. Quantum theory 

The quantization of the model is straightforward. A quantum representation of the observable algebra lfT"3")) is 
provided in the auxiliary Hilbert space H a ux = L2[SU(2) 2n , dUf], where dUf is the Haar measure. That is, states 
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have the form ip(Uf). The operators Uf are diagonal and the operators Et are given by the left invariant vector 
fields on each SU{2) element. Consistently, operators Ef-i act as right invariant vector fields. Not surprisingly, the 
operator associated to the volume V t turns out to be the standard LQG volume operator Q. 

States which are solutions of the Gaufi constraint ([bfjl are labeled by SU(2) spin networks on the graph A* . The 
dual triangulation A* is characterized by a node for each tetrahedron and a link for each face of the triangulation A„. 
In the dual triangulation only four-valent vertices appear, which intertwine spin-j/ representations Dj,(U) associated 
to the holonomies around the links If. The basis of these spin- network states is labeled by t t ), where it denotes 
the intertwiner quantum number at the given node. The spin network basis states are explicitly given by 

= {Uf\jf,H) = ®/ D^f\U f ) • ® t it , (17) 

in which "•" indicates the contraction of the indices of the D^ f >(U) matrices with the indices of the intertwiners i t . 

This constrained model can be directly quantized a la Dirac. The Hamiltonian constraint can be defined a la 
Thiemann, rewriting (|15l) in the form 

C t = ]T e fff " Tr[U ff >Upi{U r ,Vt}]ttO (18) 
ff'f'et 

and then defining the corresponding quantum operator by replacing the Poisson bracket with the commutator. Oth- 
erwise, there exists a second possibility: to implement directly as a quantum operator the constraint (|15p rescaled by 
the volume Vt 

C t = V t Ct = £ T±[Uff,E f ,E f ] « . (19) 
ff'et 

This is exactly the early proposal to perform the quantization of the Hamiltonian constraint in LQG, see [21| . Physical 
states are those annihilated by the quantum version of (I19|) . i.e. we have to impose Ctip = 0. The Hamiltonian 
constraint on the whole space £ is obtained combining the set of n Hamiltonian constraints by means of the lapse 
function N = {N t } at each node. Namely, there is C(N)ip = J2t N t^t tp = 0, ViV. 



IV. DIPOLE COSMOLOGY 



A topological three-sphere can be constructed by gluing together the boundaries of two three-balls (more details in 
Appendix |XJ. The boundary of a three-ball is a two-sphere, and these two two-spheres are identified, with opposite 
orientations. The boundaries of the two balls define then an "equator" of the three-sphere. Accordingly, a three- 
sphere can be triangulated by gluing together two tetrahedra, along all their faces. The equator of the three-sphere 
is then triangulated by four triangles /, where / = 1, 2, 3, 4. (This is a cellular complex decomposition [22j], and not 
a A-complex or a simplicial complex decomposition [2J, [24|.) The dual graph of this triangulation described above is 
formed by two nodes joined by four links 

A* — 



From now we focus on the theory defined by this triangulation. 

The unconstrained phase space of the theory defined by this triangulation has twenty four dimensions and is 
coordinatized by (Uf,Ej). At each node there are one Hamiltonian constraint C ps and three Gaufi ones Q 1 « 0. 
But it is easy to verify that the constraints of the two nodes are in fact the same, giving a total of four constraints only. 
This brings the number of degrees of freedom down to eight. The Poisson bracket algebra between the Hamiltonian 
constraints closes in this case, as the Hamiltonian constraints at each node are actually the same one because of (|12l) . 

The Hilbert space of the quantum theory is L%\SU (2) 4 / SU (2) 2 ]. The spin network states that solve the gauge 
constraint are given by states \jf,tt) = J2, J3, J4, tii^). The action of the single Hamiltonian constraint gives in 
general C\jf,it) = X)//' ^ V f'\j f i *t) • The Hamiltonian operator acts at each node. The action of the Hamiltonian 
quantum constraint operator implies that 

Cl2b'l,j2,j3,j4,tl,'-2> = Cffllll \h + ^,h+ 2^3,^4,4,4): ( 20 ) 

e,<5=±l 
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in which matrix elements Cj can be computed using recoupling theory. In terms of the wave function components, 
a more compact notation is given by 

CtWf,*) = £ c;;:;v>(j/ + f,4) ■ (2i) 

£j=0,±l 

Notice that C^* vanishes unless £/ = for only two j. Matter fields can be simply added in this picture (see @). 

In a discrete fiducial algebra element ujj associated to the triangulation was introduced, and used to compare 
the variables of the dipole model with the FRW variables. Here we introduce a novel and more useful definition for 
u)j . Consider the Plebanski two- form of the connection ui 

E 7 H = ie 1 j K uj j Kuj k = i w e aI e abc dx b A dx c , (22) 

where cj = det(w^), and let u>j be the surface integral of this two- form on the triangle / of the triangulation. Using 
the Maurer-Cartan equation ((2J , we have 

uj) = J £ J = i J e 1 j K u J A uj k = J dto 1 = <j> u 1 . (23) 

That is, the flux of the Plebanski two-form across a triangle is equal to the line integral of uj 1 along the boundary of 
the triangle. An immediate consequence of this is that for each tetrahedron t 

5>f = 0, (24) 
Jet 

where the sum is over the triangles that bound the tetrahedron t. This is because the boundary of a boundary 
vanishes. The set of su(2) vectors uj form a natural background fiducial structure for the discrete theory, analogous 

to the uj 1 fiducial connection in the continuous theory. Notice that integral of |E| = £/ on a triangle / is (twice) 
the area of / determined by the background triad. 

In the case of the dipole, the explicit properties of the uj can be found using the symmetries. The action of SU (2) 
that transforms the triangles into one another preserves the equator of the three-sphere. The equator of SU{2) is 
formed by the 7r rotations around a direction n. This set of rotations are transformed into one another by the adjoint 
action of the group on itself (which rotates n preserving the rotation angle). A discrete subgroup of this adjoint action 
sends therefore the triangles into one another. Under this adjoint action, u 1 (which is in the algebra) transforms under 
the adjoint representation, which is the fundamental representation of SO(3). That is, the four vectors ujj are rotated 
into each other by rotations (of the / index) of the same angles. Therefore they are proportional to the normals of 
a regular tetrahedron in M 3 . The cosine of the angle between two such vectors is therefore 1/3. It is also convenient 
to take uj 1 so that the norm of ujj is one: ujujj = \u>f\ 2 = 1. This characterizes entirely the wj, up to an overall 
rotation. 



V. ISOTROPY: CLOSED FRW MODEL 

We now begin the analysis of the physics and geometry of the model. We start, in this section, with the restriction 
of the model to the homogeneous and isotropic sector. This analysis was already presented in [6(, but we show here 
that it can be strongly ameliorated by using a different ansatz on the relation between the model variables and the 
scale factor. Anisotropies and inhomogeneities are discussed in the following sections. 

A. Classical framework 

Let c(t) and p(t) denote the isotropic connection variables used in LQC [illl3j. They are related to the scale factor 
a and to its time derivative d by c = (& + 1)/2 and \p\ = a 2 (see Section Ql]) . The phase space of an homogeneous and 
isotropic cosmology is two-dimensional and the basic Poisson bracket is {c,p} = I. 

In [6j, the isotropic phase space variables (c,p) are identified with a subspace of the dipole model phase space by 
the relations 



Uf — exp (coj^ti) , Ef = pto^Ti. 



(25) 
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These relations define an embedding of the isotropic geometries in the phase space of the model. Notice that this is 
different from what is done in (l5j . where a projection of the anisotropic degrees of freedom of the Bianchi I model 
to those ones of the FRW model was defined. (See also [25| for a discussion about superselection of isotropic FRW 
states of LQC within the anisotropic Bianchi I setting.) 

The second of the relations (1251) follows immediatly from ([9]), (|23|) and the identification of Ef with the flux of the 
electric field E though / (Section IIII[) . The discussion of the first relation is more delicate. Recall that Uf is the 
holonomy of the connection along a dual link e/, 

U f ^Vexpf Ai(e f ) n de<}. (26) 

Je f 

Let us assume for simplicity that the identity in the group G is identified with the center of one of the two tetrahedra. 
Then the four dual links e/ run along one-parameter subgroups of G. Now consider a point g on the SU(2) group 
manifold. Because of the flatness of the Maurer-Cartan connection uj, the holonomy of uj from the identity to the 
point xp does not depend on the particular path which has been chosen to connect the identity with g. We can write 
g = exp(a / r / ), in which a 1 = (^{g) is a vector in the internal space, and then parametrize the path 7 from the 
identity of the group to g by 

g(s) = exp(soV) , (27) 
in which s € [0, 1]. This particular path defines an abelian subgroup of SU{2), and thus the holonomy reads 

U f = Vexp I wV = g(s) . (28) 

That is: the holonomy of uj from the identity to g is precisely g. Consider now the tangent to the path 7 at g, and 
denote it v a G T g (T g is the tangent space of the group). Since we have a metric structure °q a b = ^i^i defined on 
the group manifold, we have also the corresponding one- form n a = °q a b v b . The contraction of ui and v a ds is given by 

uj{v a ds) = riae^T 1 ds . (29) 

On the other hand, by means of ([1]) and (|2"7)) . we find 

uj(v a ds) = g^dgiv* ds) = aVds . (30) 



Hence 



a 1 . (31) 



Now, take g to be the intersection between the dual link e/ and the triangle /. By symmetry, the link e/ and the 
triangle / are orthogonal. The last equation shows that a 1 (g) is nothing else than the the normal n a to the triangle 
/, with the index raised by the fiducial triad. But uif is precisely the integral of efn a on the triangle, and we expect 
it to be proportional to the value of efn a on g. Therefore we can conclude that 

U f = exp (kujjn) , (32) 

for some k. Following [11], differently from what discussed in Q we make the ansatz 

Uf = exp ((c + a) ujfTj) , (33) 

where a is a numerical constant determined by the spatial curvature, which we fix below. We refer to [ll| for a 
discussion of the rational of this identification. 
The Gaufi constraint reads 

/ / 

and is automatically satisfied because of (|24p . The scalar constraint (fTS"]) reads 

C = p 2 J2 Tr [e c ^ T 'e a ^ Tj e~ auJ f' TK e- CUJ f' TL u J s f , uf t s t m ] « 0. (35) 



8 



Using the relation e CM f Tl = cos(c/2) + 2lo^ti sin(c/2) the expression above can be simplified. Firstly notice that the 
term containing two r-matrices is zero because of the Gaufi constraint. The terms including the traces of three and 
five r-matrices are also zero because of the internal product of two equal vectors luj and of the Gaufi constraint. 
Moreover the term proportional to the trace of six r-matrices identically vanishes. (See Appendix IB"! for details about 
traces of r-matrices.) Finally, as already discussed, each identifies a normal to the center of the face / of the 
tetrahedron, the angle between them given by cos9ff> = ujuj, = 1/3. Collecting these considerations and using the 
relations 

E^H' - E cos %' = °. EW) WO = E cos20 //' = I E m 2 m 2 = e, (36) 

we obtain the final form of the scalar constraint 

C = —p 2 (cos(c - a) - 1) « 0. (37) 
6 

Such a constraint describes the dynamics of a (curved) triangulated closed FRW cosmological model. The ordinary 
classical dynamics can be recovered as soon as the limit of small connection |c| <C 1 is taken into account. The 
constraint (|3"T)) rewrites as 

17 / c 2 \ 

C= —p 2 I csina- yCOsa-2sin 2 (a/2)J + C(p 2 c 3 ) w 0. (38) 

This is the FRW classical constraint Cfbw ^ —p 2 c(c — 1) as soon as a is fixed at a suitable value a* . In fact, the 
canonical transformation c — > ac + b, p — > p/a (a, b G R) with cos a* = (9 — \/l7)/8 transforms the two constraint 
into one another. The appearance of the ordinary dynamics for small values of the connection c is in agreement with 
the claim that a coarse triangulation well approximates the classical theory for a low-curvature space-time. 

The expression (|3"T|) can be regarded as the effective Hamiltonian constraint of the standard LQC. The effective 
formulation of LQC is usually obtained by "polymerizing" the classical model: the connection is replaced by its 
exponentiated version, to reflect the fact that the connection operator does not exist in the LQC Hilbert space 
[26j |. The effective LQC dynamics can be also obtained by using the methods of geometric quantum mechanics 12711 . 
Replacing c — > sm.{vc)/v (where v denotes the polymer scale, fixed by the "minimal area gap" argument [26l . |28| ) 
in Cfrw-i one obtains the constraint Q37[) for a = a*. Here, instead, the "polymerization" is a consequence of the 
existence of the triangulation, and a conventional quantization of the truncated theory. In the quantum theory, a 
WDW difference equation arises without adding an input from LQG by hand: the lattice model leads directly to a 
difference evolution equation, without recurring to the "minimal area gap" argument. 



B. Quantum framework 

The variable c multiplies the generator of a U(l) subgroup of the compact group SU(2) 4 . Therefore it is a periodic 
variable: c € [0,4-71"]. The kinematic Hilbert space Hi SO of the theory is thus L 2 (S' 1 , dc/Aii) of square integrable 
functions on a circle. Eigenstates of p (from now on we drop the hat in order to simplify the notation) are labeled 
by an integer \i and read (c|/x) = e iAIC / 2 . Wave functions ip(c) are decomposed in a Fourier series of eigenstates of p, 
labelled by an integer \x 

V^c) = ]T ^ . (39) 

n 

The fundamental operators on this representation are p and exp(zc/2), whose action on generic states reads 

p = n/2 \n) , exp(ic/2) \ix) = |/x + 1) . (40) 
In particular, the operator sin(c/2) acts on as 

sin(c/2)|/z) = i(| M +l)-| M -l)) (41) 
2,1 

and it is then immediate to obtain the action of the composite operators sin(c/2) cos(c/2) and sin 2 (c/2). 
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By following this path, we can construct the quantum constraint operator corresponding to (1371) . In particular, in 
order for the states ip(c) to be physical states, the coefficients tpn = (c\ p) have to satisfy the recurrence relation 



where 



D+(p) (c| /i + 2) + D°(p) (c| pi) + D~(jjl) (c\ p - 2) = , 



17 

D+(p) = — p 1 ( -isin(a/2)cos(a/2) - -sin 2 (a/2) + -cos 2 (a/2) 



6 

^ 2 Q-sin 2 (a/2) 

— /z 2 f isin(a/2)cos(a/2) - isin 2 (a/2) + icos 2 (a/2) 
6 V 2 2 



(42) 
(43) 



In (|42[) we have chosen the simplest normal ordering of posing multiplicative operators on the left of derivative ones. 
Equation (|4"2")l has the structure of the LQC difference equation. 

VI. ANISOTROPY: BIANCHI IX MODEL 

In this Section we free more degrees of freedom than the sole scale factor and we describe the triangulated version 
of a homogeneous but non isotropic space-time. Since the space topology is S 3 , we are therefore dealing with the 
Bianchi IX cosmological model. Bianchi IX (or the Mixmaster Universe [19|]) is the most general homogeneous model 
with this topology: its physical importance relies in describing a generic (classical) solution of the Einstein equations 
toward a space- like singularity via the BKL scenario Q- 



Classical framework 



The dynamics of an anisotropic (homogeneous) cosmological model is described by three scale factors a/ = aj(t), 
which identify three independent directions (in the time evolution) of the Cauchy surfaces. In the connection for- 
malism, relaxing the isotropy condition corresponds to consider three different connections cj = ci(t) and momenta 
pi = pi(t). Our model can be then extended to an anisotropic setting by demanding that the variables of the theory 
are given by 



Uf = exp (c 7 o;^r/) exp (a W/Tj) , Ef 



p T UfTl. 



(44) 



Notice that there is a main difference between the holonomies (|44j) and those used in LQC in the anisotropic context 
[Til [T^|. In the standard formulation of an anisotropic LQC model the basic holonomies hj, are directional objects, 
computed along edges parallel to the three axis individuated by the anisotropies. They read 



hj = cxp(c (/) w (/ )) = cos(<//2) + 2u>i sin(c / /2). 



(45) 



On the other hand, the variables Uf are non-directional objects, because the four faces of the triangolation do not 
have any special orientation with respect to the three isotropy axes. The connection components are summed over 
and they are thus independent on the /-direction. The Uf are in fact group element of SU(2) that depend on the 
face /. Explicitly, they are given by 



exp (c WfTi) — cos 



i, J 



C UJ 



+ 2 TTTT sin 



cWf 



cos Pf + 2p I fTi sinp/, 



where denotes the unit vector p^ 



Pf/pf defined by 



Pf 



pf 



= jY,(p I f) 2 = l\c I » I f\- 



(46) 



(47) 



Let us now analyse the dynamics of the anisotropic model. As above, the Gaufi constraint does not carry out any 



information since it vanishes because of the construction of the vectors uj j. Explicitly, it is given by 



>,.,(/> = 



(48) 
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The scalar constraint of the triangulated Bianchi IX model reads 



C = J^Tr [(cos/?/ + 2/3/T/sin/?/) (cos(a/2) + 2w/Tj sin(a/2)) • 
//' 

• (cos(a/2) - 2uf,T K sin(a/2)) (cosp// - 2^,r L sin t s t m ] 

where the £j are defined by 



0, (49) 



(50) 



The difference with respect to the usual LQC constraint is evident. Let us note that, as a result of taking Uf as 
defined above, no terms drop here in virtue of the Gaufi constraint, as they did in the isotropic case. In fact, terms of 
the kind gfgfUifOJf (where <?/ are generic functions) appear and their sum in /, /' no longer vanishes. It is convenient 
to define 



Xff 



ujpj, 



Pff'=Uftf>> Iff 



a ff 



S ff = PfPj' 



and 



A = / 5A£A£, 
D=pAu}A£, 



B = lu A lo A £, 
E = p A p A u), 



C = u A cj A p, 
F = ujA£A£. 



(51) 
(52) 



where the convention is A//// = pf A£/ A £/' = euKpi^f^f and so on. The scalar constraint is given, in terms of all 
these quantities, by 



C = < sin p fi cos p f 
1 



A ff'f cos («/ 2 ) + sin(a/2) cos(a/2) Pfylf'f 



'flff + \xpf°ff + Pfflff 



Pfflf'f ~ oXff'Vff + 2 sin (a/2) - lff B ff , f - - lff ,B ff , f , + a fr C frr + cos%M///,/ 



+ COS Pf COS pf 

+ sin p f cos p f 
1 

+ oXffVff 



COS tiff Off 



-Off cos 2 (a/2) 



2sin(a/2)cos(a/2)F / , / , / +sin 2 (a/2) [PtfPtr - PtfPff 

\ A ff'f cos 2 (a/2) + sin(a/2) cos(a/2) (jffPff ~ IffPf'p ~ Iff'Pff + VfPff ~ \xff>°ff>+ 

cos 2 (a/2> 



*iu J l"/2) ( \fiffDfff - ^PffDfff+OffCfff + X ffFff f 



smpfsmpf 



1 



1 



1 



If f'lff ~ Ifflff ~ o 5 //" 7 //' + sin(a/2) cos(a/2) -jffDfff - -rJffDfff - 2o f fE f ff + Xff' A f'f'f+ 



+ ^lff D fff ~ llff D fff ~ Xff A fff j + 2 sin ( a / 2 ) \^2 E fff F ff'f ~ 2 Cfff ' Aff ' f + aff ' ( cosd ff' S ff' + 

-Xff Xff) + Xff (iff'Pff ~ IffPff) + Xff (iffPff ~ IffPf'f) ~ ^XffXf'fVff ) } ~ . (53) 

The dynamics of the Bianchi IX model, in the triangulation formalism, is then expressed in terms of geometrical 
quantities as angle-like (|5"Tj) and volume-like (|5"2j) quantities. As can be checked by direct calculation, the isotropic 
constraint (|37|) is recovered as soon as c% = c 2 = c 3 and pi = p\ = f?3- In this case all the terms in (|52"j) vanish because 
the internal wedge product of two equal elements Uf. Furthermore the terms containing any single angle (|51|) are also 
zero because of the Gaufi constraint. 

As in the isotropic case, the constraint (|53l) describes a discrete (triangulated) dynamics. Such a dynamics can be 
regarded as an effective one in which all the outcomes of the triangulation are abridged. The standard Bianchi IX 
model is recovered in the limit of small connections, i.e. as soon as |c/| <C 1. In this limit, which corresponds to take 
\pf\ <C 1, the expression (|53[) is rewritten as 



C = p*p J (const + c K + c K c L ) + 0(p 2 c 3 ^ 



0. 



(54) 



In order to make the notation the easiest as possible, constants (depending also on a) are not written in (I54[) and 
the contraction of the free indices is understood. The isotropic constraint is recovered if a basic canonical 
transformation is performed, i.e. cj — > c/ + dj and p 1 — > p 1 + b 1 for d 1 \b T € R 3 . 
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Summarizing, the constraint (|53| describes the effective LQG dynamics of an anisotropic homogeneous space-time 
having the spatial topology of a three-sphere, i.e. the Bianchi IX cosmological model. In the triangulation formalism 
the anisotropics (different scale factors or different connections) are described in terms of angles (|5Tj) and (internal) 
volumes (f52"j) characterizing a generic dynamical tetrahedron. The contact with the standard (classical) dynamics is 
recovered as soon as the flat space-time limit \ci\ -C 1 is taken into account. 



B. Quantum framework 



The fact that the dipole model includes Bianchi IX suggests that we may use it in order to construct a loop 
quantization of Bianchi IX. We analyze this intriguing possibility in this section. There are two possible ways for 
doing so. One is to start from the (c 1 ,pi) variables, and define the dynamics with the Hamiltonian operator (|53|) . 
where all quantities are taken as functions of (c 7 ,pi) via equations (|4"T)) and (|50[) . The second possibility is to use 
the full quantum theory of the dipole cosmology model, defined in Section IV Bl Neither of these approaches gives a 
simple quantum model, as far as we can see, since the resulting dynamical equation is complicated in both cases. 

Let us begin with a basis of states |c 7 ) that diagonalize the variables c 1 , and write the quantum state as ipic 1 ) in 
this basis. The p 1 operator is then —id/dc 1 . For each face /, consider the three scaled variables (l47l) and the three 
scaled momenta (|50[) . The Hamiltonian operator (1531) is entirely written in terms of the four operators sin pf, 
cos pf and p^/pf. These operators are easily written in this representation. The first is 

£/ = ' ( 55 ) 

while the others are diagonal. Thus, (|53l) yields a Hamiltonian constraint on this Hilbert space, and defines a 
quantization of the Bianchi IX cosmological model. 

Does this quantization yield discreteness? Fix for a moment one of the four faces, say /. Observe that p^ enters the 
dynamical theory only as an argument of the group element Uf in (144|) via Uf = exp(2/5yT/) exp(aw^r/). Therefore 
the physically relevant domain of pj appears to be given by 

\Pf\<^ (56) 

that is, adding Att to \pj\ gives back the same configuration of the Uf variable (for a fixed /). Accordingly, the 
configuration space is compact, and we expect the conjugate variable pi to have discrete spectrum. This is reflected 
in the fact that the norm of Pf enters the Hamiltonian only via sin pf (or cosp/). It seems therefore reasonable to 
expect that what is going on is something analogous to what happens in standard LQC, namely we can chose a scalar 
product where momenta have discrete spectrum, the operator sin pf is well defined while pf is not, and the dynamics 
is restricted to a discrete set of eigenstates. These observations are intriguing, but incomplete, because there are four 
distinct faces /, and four distinct group elements Uf, and generically, there is no shift on the c 1 that gives back the 
same value for all four group elements. In other words, the map c 1 —> Uf defined in (|44[) maps R 3 into the compact 
space [SU(2)] 4 , but its image is not necessarily compact. We leave the development of these reasonings for future 
investigations. 

The second alternative is to use the dipole model Hilbert space "Hd, formed by states ip(Uf), and spanned by the 
discrete \jf,i) basis. The reduction of a quantum theory to a sector where some degrees of freedom are frozen (here 
the inhomogeneous ones) is, in general, a rather non trivial step Q. If we have given the Bianchi IX Hilbert space 
Wix formed by functions ipic 1 ) and the c 1 —> Uf map (I44|) . then a natural projection n : Hd — > Hix is simply defined 
by i/j(cj) — ip(U /(c 7 )). But in order to map v : Hix T~ld we need to choose a way to compute c 1 out of Uf. One 
natural possibility is for instance to define 



c 



p (I) 

'=E £777- («> 



Using this, we can map any state in Hzx to 7-Ld- For instance, we can expand an eigenstate exp(i/i / c / ) = Yli (c 7 
of p 1 on the spin network basis, considering Uf = Uf(pf) and formula (1171) . via 

n07.il/) = / du f \®f Djf (Uf)-i-i] expfiXXE \A . (58) 
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Using these coefficients, the Bianchi IX Hamiltonian operator can be simply defined by Cjx = ttCv, where C is 
the dipole Hamiltonian constraint operator, defined in the Section IVIII This operator will be studied more in detail 
elsewhere. 

VII. INHOMOGENEITY: PERTURBATION OF BIANCHI IX 

Let us finally come to the full geometrical interpretation of the "dipole model" . The remaining large-scale grav- 
itational degrees of freedom captured by the dipole dynamics are necessarily inhomogeneous. The easiest way to 
characterize them is to consider an expansion of the gravitational fields in tensor harmonics, and identify them with 
the lowest order term of this expansion. In other words, we can imagine that we want to describe only some low order 
term of a mode expansion of the geometry over the three sphere, and that the dipole variables read out these terms. 
In [2t|, Regge and Hu considered a similar mode expansion around Bianchi IX, using Wigner D-functions. (See also 
[3d||7) Here we follow their approach, adapting it to the first order formalism we are using. 

We begin by recalling the Regge-Hu expansion [2j|. The Wigner D-functions D 3 a , a (g{x)) determine a basis of 
functions of the symmetry-group of the model. Recall that we can use group elements g(x) to coordinatize the physical 
space that has the S 3 topology. Let q° b (x,t) be the metric of an homogeneous space. A general perturbations to the 
three-metric h a b(x, t) = q a b(x, t) — q^ b {x, t) can first be translated into a matrix of space-scalars hu(x, t) by projecting 
it on the invariant one-forms, that is 

h a b(x, t) = hu(x, t) Uj(s) <i){{x); (59) 

Regge and Hu decompose these scalars in terms of definite angular-momentum components of the three-metric 
h J jj(x,t), labelled by spin and magnetic numbers {j, a}. 

h IJ (x,t) = J2 Kj(x,t). (60) 
These can be expressed in terms of Wigner D- functions D 3 , (g(x)) (whose form is recalled in Appendix via (see 

h%(x,t)= j2 rtr'WB'a'M*))-, ( 61 ) 

a'=—j 

The time dependent amplitudes h 3 j" a (t) represent, at fixed j, ( 2 j + l) 2 inhomogeneous degrees of freedom. They are 
governed by a set of coupled differential equations, studied in [3l| . 

At a first sight, the dependence of the scalar harmonic function h 3 j"(x,t) on the modes j, a in equations (|60H61j) 
seems to be incompatible with the Peter- Weyl decomposition [32| of hu(x,t). However, the choice of not summing 
over the j, a labels contracted with the Wigner D functions relies on the evidence, due to Einstein equations, that 
j, a states can be de-coupled from a' states. This leaves as a possible choice to fix perturbations of the metric with 
definite j, a. Indeed, only the a' states are mixed by the action of the derivative operators and thus by the linearized 
Einstein tensor, as explained in Appendix [Cj This can be understood moving from the action (|C9[) of the invariant 
operators on the Wigner D functions and expressing it in terms of derivatives vector fields on the Euclidean space E A 
in which S 3 can be embedded. 

Let us adapt this formalism to the first order variables we use. We start with the triads. Let us write a generic 
perturbed triad Ef(x, t) as the sum of the background triad e/ field and a perturbation. 

Ef(x,t) =e a I {x)+rAx,t). (62) 

It is convenient to project the perturbation on the background triad 

r i (x) = ^u(x,t)e a J (x). (63) 

Following Regge and Hu, we write this as a sum of components of definite j and a quantum numbers 

i; I j{x,t)=Y J ^ I j{x,t) , (64) 

j a 
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where 

4j{*,t)= E 47\t)DlMx)). (65) 
The same can be done for the connection 

"> *) = + *) ■ ( 66 ) 

Project </3^ on the invariant one-forms 

^(z)=^ J (M)KfaO. (67) 
Expanding this in components of definite j and a quantum numbers gives 

=E</^(M) , (68) 

j a 

where 



^(x,t)= £ ^L{t)Di, a {g(x)). (69) 



The (Vjaa'W'^lia'W) are tne time-dependent expansion coefficients that capture the inhomogeneous degrees of 
freedom. They, are given by matrices in the internal indices /, J, labeled by the spin j that runs from j = 1/2 to all 
the semi-integers numbers, and the corresponding magnetic number a. 



A. Dipole model 

Suppose we restrict the geometry by assuming that the matrices (<pja a i{t)j V'/aa'W) are diagonal in the internal 
indices /, J and it is different from zero only for lowest nontrivial integer spin j = 1 and for, say, a = 0. That is, we 
restrict to the components 

where a = —1,0,1. Then we have that the inhomogeneities are determined by precisely nine plus nine variables 
fait), tp a (t), namely nine degrees of freedom. Assuming that the Gaufi constraint reduces the degrees of freedom by 
three (see below) , we obtain six degrees of freedom, which is the number of degrees of freedom captured by the dipole 
variables. Thus, such a geometry is entirely captured by the six degrees of freedom of the dipole model. Therefore, 
we can interpret the six extra degrees of freedom of the dipole model (beyond anisotropics) , as a description of the 
diagonal part of the lowest integer mode of the inhomogeneities 

If we do so, we can relate the variables of the dipole model to the quantities (f a (t) , tjj a (t)) . Notice that the one 
forms Cj 1 no longer satisfy the Maurer-Cartan equation ([2]). The fiducial algebra-elements ujj are therefore perturbed 
as well. At first order, for a generic perturbation, let us define 

= \ [ e 1 jk & A Q K = coj + J e 1 JK lo j A tp K — ^ + £ <P$Z, («) (71) 

* jaa f 

where 

^Tl = j e 1 JK Di a ,u^u L . (72) 
In particular, if we restrict to the diagonal j = 1, a = case, 

where 

4>), a = J e 1 j K Dl a u J Au K (74) 
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are fixed coefficients. Then the relation with the dipole variables can be written as 

Ufic 1 , tpi) = exp (duifTi) exp [aui^n) , (75) 
which replaces (|44l) . Similarly, we can write 



Ef = [ (c? + ^i)e ab cdx b A dx c = E) + 2 f $ u e J KL u K A i> L = Ej + 2 Y] ifitfj aa , 
Jf Jf 



(76) 



In particular, if we restrict to the diagonal j = 1, a = case, 

^ = Ej + 2^ 4 a . (77) 

The Gaufi constraint no longer identically vanishes. It can be split into two parts: the homogeneous and the 
inhomogeneous terms 

g i = J2p^; ) + 2 ^ Etf.«« - ( 78 ) 



The first part is the constraint which appears in (|48|) within the Bianchi IX framework and vanishes identically because 
of the Stokes theorem. The second gives three conditions on the inhomogeneous perturbations to the electric fields 
Ef- 



VIII. CONCLUSIONS AND PERSPECTIVES 



We have studied a triangulated (loop) quantum cosmology model, and analyzed the geometry it describes. In 
particular, the model of a dipole S'f/(2)-lattice theory, triangulating a topological three-sphere by means of two 
tetrahedra, has been related to a Bianchi IX cosmological model perturbed by six inhomogeneous degrees of freedom. 

This work is mainly based on the case of the dipole cosmology, for which the algebra of Hamiltonian constraint 
closes. Higher modes of the mode expansion can be captured with finer triangulations, but these require a more work 
for a consistent definition. 

The truncation we use relies on spherical topology. It is likely that the technique we have used could be extended 
to general compact spaces, but we have not studied this issue. For the extension to non-compact spaces, on the 
other hand, the problem is delicate. To use the same technique, one should partition a non-compact spatial slice with 
"fiducial boxes", as is done in LQC, but this would reduce the inhomogeneities that the model can capture. Doing 
so, on the other hand, would be interesting for the comparison with LQC. 

We have only discussed the dynamics generated by the Euclidean part of the Hamiltonian constraint and with a 
fixed value of the Barbero-Immirzi parameter. The full Hamiltonian constraint, the possibility of using of a generic 
value of the Barbero-Immirzi parameter j3, and the Lorentzian theory, will be considered elsewhere. 

Finally, notice that the results of this paper, together with the link between LQC and spin- foams derived in [4lj |. 
might provide a path to connect cosmological spin-foam models from the cosmological sector of LQG. 
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Appendix A: The three-sphere 



We collect in this Appendix some useful information on the structure of a three sphere, and a few other useful 
formula. In IA II we illustrate the topology of the three-sphere and its relation to the Hopf fibration. In subsection 
IA 21 using the isomorphim of groups S 3 /{1, —1} ~ 50(3), in which the subgroup {1, —1} is the kernel of the group 
homomorphism h : S 3 —> S0(3), and identifying with SU(2) the double cover of 5*0(3), namely Spin(3), we clarify 
which is the role of the group SU (2) in describing the topology of S 3 by means of the Hopf fibration. In subsection 
IA 3[ we describe how the SU(2) symmetry structure enter the definition of the Cartan one-forms. 
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1. Topology of the three-sphere: Hopf fibration 

We illustrate here a concrete realisation of the Hopf map [33] acting on S 3 , which defines a fibration with fiber space 
S 1 and the base space S 2 . In what follows we identify R 4 with C 2 and R 3 with Cxi. Taking (x\, x 2l x 3 , X4) € R 4 
and (^1,^27^3) G R 3 the identification is achieved respectively by means of zq = x\ + ix 2 and z\ = x 3 + ix 4 for R 4 , 
and z = £1 + it; 2 and y = £3 for R 3 . A unit three-sphere embedded in R 4 , as (xi) 2 + (x 2 ) 2 + (x 3 ) 2 + (a^) 2 = 1, reads 
in terms of the complex variables |zo| 2 + \zi\ 2 — 1- A unit two-sphere is identified with the subset of C x R such that 
\z\ 2 + y 2 = 1. Let us parametrize S 2 by means of the expression (£i) 2 + (£2) 2 + (£3) 2 = 1 an d take a natural bundle 
projection by the Hopf map it : S 3 — > S 2 

£i = 2(x x x 3 + x 2 x 4 ), (Al) 

£2 = 2(x 2 x 3 - x x x 4 ), 

£3 = -{xi) 2 -{x 2 ) 2 + {x 3 ) 2 + {x4) 2 . 

With this mapping the relation (£ x ) 2 + (£ 2 ) 2 + (£ 3 ) 2 = (O^i) 2 + (x 2 ) 2 + (x 3 ) 2 + (xi) 2 ) 2 = 1 holds. 

We can now consider a stereographic projection [§3, HH of a point in the southern hemisphere S 2 . C S 2 from its 
north pole and label its coordinates as U, V. On the complex plane C containing the equator of S 2 C C x R, we 
consider the function T = U + iV, which is within the circle of unit radius on the complex plane and is given by 

j, _ + i£ 2 = x x + i x 2 = zo ^ A2 ^ 
1 — £ 3 x 3 + i X4 zi ' 

with G S* 2 . For a A G U(l) such that |A| = 1, we can see that T is invariant under (z , ^i) —> (Azq, which 
are both points in S 3 . Since the set of complex numbers A with |A| 2 = 1 form the unit circle in the complex plane, 
it follows that for each point m € S 2 , the inverse image 7r _1 (m) is a circle, i.e. 7t~ 1 (to) ~ 5 1 . Thus the 3-sphere is 
realized as a disjoint union of these circular fibers. 

Similarly, the stereographic coordinates P, Q of the northern hemisphere projected from the south pole read 

s = e-if = x 3 +i Xi = z 1 

1 + ^ Xi +1X2 Zo 

with f € S%. On the equator, S^nS!, one finds T = S 1-1 . 

The fiber bundle structure is then given by the following local trivialization 1 : on the south hemisphere one defines 
: 7r _1 (S' 2 ) — > S" 2 x U(l) as given by (z , z{) — > (z /zi, Zi/\zi\) while on the north hemisphere the quantity 
(p^ 1 : tt^ 1 (S^_) — > x U(l) is given by (zq, Z\) — > {z\/zq, zq/\zq\). This construction clarifies that S 2 is the base 
space of S 3 . 



2. 5(7(2) description of the three-sphere 



A further geometric interpretation of the Hopf fibration can be obtained considering rotations of the two-sphere in 
a three-dimensional space [35|]. We recall that Spin(3) is the double cover of the rotation group SO(3) and that is 
diffcomorphic to S 3 . The spin group acts in a transitive way on S 2 by rotations. The stabilizer subgroup 2 of a point 
of S* 2 is isomorphic to the circle group, from which it follows that S 3 is a principle circle bundle over S 2 . This is 
exactly the Hopf fibration. 

Concretely, this can be seen identifying Spin(3) with the group of unit quaternions Sp(l): a point (xx, x 2 , x 3 , 2:4) € 
R 4 is interpreted [36] as a quaternion q € H by q — x\ + \ x 2 + j x 3 + kx4 and the three-sphere is then identified with 
the quaternions of unit norm, namely qq = \q\ 2 = (xi) 2 + (x 2 ) 2 + (x 3 ) 2 + (X4) 2 = 1. Now a vector p in R 3 can be 
identified with the imaginary part of a quaternion, i.e. p = ipi +jp 2 + hp 3 , and the mapping p — > qpq is a rotation in 



1 Each local trivialization is well defined on each chart, while on the equator they become 0_ 1 : (zo, Zi) — > (zo/zi, \/2zi) and (f)^ 1 : 
(zo, zi) — > (zi/zq, \/2zo): the transition function on the equator is t_|._ (C) = (v / 2zo)/(v / 22 1 ) = e E/(l). While following a path 
around the equator, the transition function t_|-_(^) crosses the unit circle in the complex plane once. Thus the (7(1) bundle S 3 % S 2 is 
characterized by the homotopy class 1 of 7ri((7(l)) = TL. The Hopf fibration is therefore locally a product space, but it is not a trivial 
fiber bundle, as S 3 is not (globally) a product of S 2 and S 1 . 

2 For every element x of a set X, the stabilizer subgroup of x is the set of all elements in G that fix x. 
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R 3 , and thus an isometry as \qpq\ 2 = \p\ 2 . Morover it is not hard to check that this mapping, as a rotation, preserves 
orientation. The q are then unit quaternions provided by the group of rotation in R 3 , with opposite elements — q 
undergoing the same transformations. The set of unit quaternions q which fixes a unit imaginary quaternion have 
the form q = ( + vp, with (,i/£R such that |£| 2 + W\ 2 = 1 {i-e. a circle subgroup). Taking then p = k, one can 
define the Hopf fibration via the map tt : q — > gkq. The image set by the Hopf map is made of (still unit) quaternions 
provided only by imaginary parts: these points lie on the two-sphere S 2 . As a quaternion q = x\ + i X2 + j £3 + k £4 
can be recast in terms of a 2 x 2 matrix, 



one can identify the group of unit quaternions with SU(2) and represent imaginary quaternions by the skew-hermitian 
2x2 matrices, which are isomorphic to CxK. The fiber for a given point of S 2 consists now of all those unit quaternions 
whose image via the Hopf map n : q — > qkq aims there. These points are easily recognised to belong to a circle. A 
direct way of seeing it is to consider that the multiplication by unit quaternions is equivalent to realise a composition 
of rotations in R 3 : for instance multiplying by q t — e kt corresponds to a rotation by 2i around the z axis — varying 
t one draws a great circle of S 3 . As long as a generic base point /3 — (l,m,n) is not the antipode (0,0,-1), any 
quaternion qp (associated to the base point /3 by means of ir~ 1 ) not having component on the z axis, say for instance 
qp = (1 + n, —m, I, 0)/-\/2(l + n), will continue aim there after multiplication by q t : the fiber of the base point /3 is 
therefore given by quaternions of the form qg - qt, which does not represent merely a topological circle, but a geometric 
one [33| • For the base point (0, 0, —1), whose associated quaternion can be taken to be the i axis, the fiber is simply 
(0, cost, — sint, 0). 



The unit quaternions, once expressed by means of the unimodular (because of the qq = 1 condition identifying a 
S 3 ) matrices (|A4|) . can be in turn represented by means of Euler angles, which are another set of coordinates {9, cf>, ip} 
on the three-sphere such that 9 and cf>, as usual, are in the range 9 £ [0,ir), (j> £ [0, 2tt), while for the Euler angle ip 
one must set i/j £ [0,47r) in order to achieve a simply connected covering 3 space with the topology of S 3 . By these 
one q assumes the form 



In terms of Euler angles, isometries of the three-sphere are well understood as left (and right) translation generated by 
SU(2) representations. Recall that a left (right) translation is a mapping \i : S 3 — > S 3 which take a point P in a point 
point P', where q(P) — /iq(P') (and respectively q{P) = q(P') £«)■ These are really translations, as the only mapping 
having fixed points is the identity. Otherwise, a mapping having fixed points can be defined by q(P') = /i(P)/i _1 . 
This latter is clearly a rotation. 

The right and the left translation can be used to define differential operators on S 3 . For instance a right translation 
can be represented by a hermitian matrix v and a parameter v, for which g(P') = q(P)e lvv or by a differential 
operator K and a parameter v, such that g(P') = [e lvK q](V) . The compatibility condition for those two representation 
is q v = K-q. 

Similarly, a left translation is defined by v q = Kq, where using (|A4|) and previous definition of left and right 
translation, it is possible to see that K{x\, X2, £3, £4) = IC(— Xi, —X2, —x%, X4). Similarly, in term of Euler angles, 
one must send (0, 9, tji) — > (—ip, —9, —</>). Right and left translation operators can then be determined from these 
definitions. If we relabel a point in R with {x\, X2, X3, x±) = {xi, X4), we can recognize, respectively, the right and 
the left translation operators to be [38( 




(A4) 



3. Differential structure of the three-sphere 



q(9, (j), V) := e" 0Tl e~ 8T2 e"^ 3 . 



(A5) 




(A6) 



3 This chart breaks down at the poles 9 = <j> = ip = 0, corresponding to the identity 1, and at 8 = <j> = and ip = 2n, corresponding to 
-1. 
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Notice that these generators 4 fulfill two su(2) algebras 

[Ji, Jj] = ie ijk Jk , [Ji, Jj] = itijkJk , [Ji, Jk] =0. (A7) 

These four-dimensional operators are all tangent to the surfaces of constant modulus, represented by {x{] 2 + ix^) 2 + 
(X3) 2 + (X4) 2 = e*, with t € M. Thus, focusing on our case t = and equating relation (|A4I) . which is written in 
cartesian coordinates, to relation (|A5[) expressed in terms of Euler angles, one finds change of variables and then 
translation operators in terms of Euler angles. 

Making a further transformation, which defines Hopf coordinates in terms of Euler angles, 

= £+ = ^ + 0, i- = \{^-4>), (A8) 

one finds the transformation from cartesian to Hopf coordinates 

x\ = sin(0/2) sin£_ , x-2 — sin(0/2) cos£_ , 23 = cos(0/2) sin£ + , 24 = cos(0/2) cos£+ . (A9) 

Transformations (IA8|) and (|A9I) allow us to think at the topology of S 3 in terms of two-torus coordinates, as both 
£i ; £2, G [0, 27r). In fact, we can consider the third component of the left and right invariant vector fields J3 and J3, 
and their combinations J3 — J3 and J3 + J3 . Their duals one-forms are now linear in c££_ and d£ + , as 

uj 3 -Cj 3 = -2 cos 2 77 <f£+, (A10) 
w 3 +tl) 3 = -2sin 2 ry^_. 



The metric on the unitary 3-sphere turns out to be expressed by 

^ 2 = 4((d77) 2 +sin 2 ry(^-) 2 +cos 2 77 (^ + ) 2 ) , (All) 

which makes clear that the trajectories of these two-parameter subgroup are constant rj surfaces and that they shrink 
in the £_ direction near the degenerate trajectory 77 — 0, and in the £ + direction near the other degenerate direction 
rj = 7r/2. In terms of the Hopf fibration, each torus, which topologically is in stead the product of two circles, can be 
thought as the stereographic projection of the inverse image of a circle of latitude of the two-sphere. 



Appendix B: Traces of r-matrices 



We give here results of straightforward calculations of traces of products of r-matrices. We remind that a basis for 
the su(2) algebra elements in the fundamental representation is given by tj = cr//(2z), in which 07 are intended to 
be the Pauli matrices. It follows that: 



Tr[r/Tj] 
Tr [t/T/Tk] 
Tr [tjTj tx tl ] 

Tr[TlTjT K T L Ts] 
Tr[riTjT K T L T S T M ] 



4 5 



— -T^KS^IJL + -^>KLtIJS + SlS^IJK + Sij€lSK 



-^S kl (5isSjm - SimSjs) + -^Su(5 ks Slai - SkmSls) - 2§ij8kl5sm 



(Bl) 



4 Rotations about the fixed point x\ = X2 = £3 = and £4 = 1 are generated by the operators = Jk ~ Jk- 
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Appendix C: Wigner D-functions 

Finally, we recall here an explicit formula for the Wigner D-functions. These are obtained from their very defini- 
tion in terms of matrix elements of the rotation operator lZ((p,9,ip) = expj— «0j x ) exp(— i0j y ) exp(— hpj z ) , in which 
jx 7 jyi jz are generators of the su(2) Lie-algebra, and can be expressed as [39( 

D> a , a = (j a'\K(<f>, 9, ^)| j a) = d? a , a {6) e'^ , 

= [CJ + - ct)l(j + a)l(j - a)!]V» £ s u+a -s fi $--Zw-«>-.» ( cos f f '^'^ (sin § f^Cl) 

in which the sum over s is over such values that the factorials are nonnegative. 

For j = I integers, the Wigner D-functions are simply related to the spherical harmonics Y^ a (9, <p) by the relation 

D l a , a {g) = (-l)-«7M(2Z + l)n-«' <t>) e M * ■ (C2) 

More in particular, the Wigner D-functions, which are the representation functions of the SO (3) and SU(2) 
simmetry-groups, can be obtained by imposing that they satisfy the following differential equations, expressed in 
terms of the Euler angles 



j _ 



d 2 
<w 2 



™t 6f e + ( £ - 2 costf^ + D^ = j(3 + 1) ^ 



^3 K, a = D J Q , a = a' Di >a , L 2 D^/ a = D J Q , Q = a & a , a . (C3) 

In the above equation we have introduced two bases of generators for the su(2) algebra and their common Casimir op- 
erator C, which is an invariant of the group. The angular momentum operators {L\, L 2 , L 3 } of the three-dimensional 
rotation group in quantum mechanics, which are the intrinsic angular momentum operators of a rigid body, are related 
to the left-invariant vector fields ef (which generate right-transformations) via the relations 

L^ielda, L 2 = ie a 2 d a , L 3 = ie a 3 d a . (C4) 

The spatial angular momentum operators {L x , L y , L z } are in turns related to the right-invariant Killing vectors 
(which generates left-transformations) by the formulae 

L x = -igd a , L y = -i$d a , L z = ~iQd a . (C5) 

In the chart provided by the Euler angles, it is straightforward to find the expression for the left invarint frames ej 
as well as that one for the right invariant vector fields £/ (which Lie drag the the frame and coframe introduced in 
Section [II]). For the frames 

d sin ib d cos 9 sin ip d 

ei = —cosyj 



e 2 = sin^;- 



89 sin 9 d<p sin9 dip ' 
d cos ip d cos 9 cost/j d 



d9 cos 9 d<p sin 9 dip 
^ = ^ (06) 

which fulfill the relation [ej, ej] = —ejj K ex- 
it turns out that the Killing vectors are in the following relation to the frames 

. , d cos^cos^ d sin 6 d 

ii = cos <p — — — + t— — = 

o9 sm9 o<p sm& dip 

= — (cos 4> cos ip — cos 9 sin <p sin ip) e\ + (cos <p sin ip + cos 9 sin <p cos ip) e 2 + (sin <p sin 9) e 3 , 

^ d ^cos#cos0 d coscp d 

89 sin 9 d(p simp dip 

= — (cos ip sin <f> + cos 9 cos <f> sin ip) e\ + (sin (p sin ip — cos 9 cos <p cos ip) &i — (cos (p sin 9) e 3 , 

d 

£ 3 = — — — = sin 9 sin ip e\ + sin 9 cosipe 2 — cos 9 e 3 , (C7) 
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from which the Lie brackets = € ij K £.k easily follow. Moreover 

[ti,ej]=£i>ej=£i<ej = 0, (C8) 

which means that ej are left- invariant vector fields (i.e. they are invariant under the left action ">" of the Killing 
vectors) and, in turns, that £j are right-invariant vector fields (i.e. they are invariant under the right action "<" of 
the frames). 

We want now to show, following [40j], that the action of the derivatives, expressed in terms of the action of the 
invariant operators, mixes only the a' states. That is what we briefly mentioned in Section I VIII while we were 
expanding inhomogeneous perturbations (intended as scalar harmonic functions on S 3 ) as a linear combination of 
Wigner D-functions. 

From the previous equations it follows that 



L+K'a = [Li+iL 2 ) Di, a = iy/(j + a>)(j -a' + l) Df a ,_ 1)a , 



L-K>« = [L 1 - l L 2 ) D> a , a =iy/(j + a> + l)(j-a>)Dl a , +1)a , 

LsDL = u>Di, a . (C9) 

The invariant operators can be rewritten in terms of Cartesian coordinates x A = {x\, X2, £3, £4} in the Euclidean 
space E 4 in which S 3 is embedded (see Appendix In stead of the three co- frames w 1 in the Euler angles chart, in 
the Euclidean space the invariant basis is given by the four forms a A , which are related to oj 1 by the transformation 
matrices Sia(xa) via 

oj 1 = 2S IA (x A )dx A . (CIO) 

Conversely, the coordinates differentials of E A are expressible in terms of lo 1 via dx A = 1/2 Si a 
The coordinates derivatives, which are the vector fields of -E 4 , are then expressed by 

2Si A (x A )ei. (Cll) 



dx A 



It turns out to be much easier, using the homogeneity of the S 3 spatial slices, evaluate the Cartesian derivatives at 
the pole xp = (x^ — 1, xi = xi = x% = 0), where the transformations matrices reduce to S[ A = —5ia thus yielding 
the relations 



d 



dx A 



-2e/. (C12) 



Equation (|C12[) allows us, by means of (|C4|) . to express derivatives in terms of invariant operators. On other hand, 
equations (|C9I) specify the actions of the invariant operators on the Wigner Z?-functions. It follows, once the Einstein 
equations for the Bianchi IX model have been rewritten in terms of Cartesian coordinates x A |29| . that the pertur- 
bations to homogeneity possess states with definite j, a and that only perturbations labelled by a' states are mixed. 
This result has been used in Section IVIII in order to de-couple j,a modes from a' modes. 
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